Sequence Alignment-Based Prediction of Myosin 7A: Structural Implications and Protein Interactions.

Myosin, a superfamily of motor proteins, obtain the energy they require for movement from ATP hydrolysis to perform various functions by binding to actin filaments. Extensive studies have clarified the diverse functions performed by the different isoforms of myosin. However, the unavailability of resolved structures has made it difficult to understand the way in which their mechanochemical cycle and structural diversity give rise to distinct functional properties. With this study, we seek to further our understanding of the structural organization of the myosin 7A motor domain by modeling the tertiary structure of myosin 7A based on its primary sequence. Multiple sequence alignment and a comparison of the models of different myosin isoforms and myosin 7A not only enabled us to identify highly conserved nucleotide binding sites but also to predict actin binding sites. In addition, the actomyosin-7A complex was predicted from the protein-protein interaction model, from which the core interface sites of actin and the myosin 7A motor domain were defined. Finally, sequence alignment and the comparison of models were used to suggest the possibility of a pliant region existing between the converter domain and lever arm of myosin 7A. The results of this study provide insights into the structure of myosin 7A that could serve as a framework for higher resolution studies in future.

Analysis of the Conformational Landscape of the N-Domains of the AAA ATPase p97: Disentangling the Continuous Conformational Variability in Partially Symmetrical Complexes.

Single-particle cryo-electron microscopy (cryo-EM) has been shown to be effective in defining the structure of macromolecules, including protein complexes. Complexes adopt different conformations and compositions to perform their biological functions. In cryo-EM, the protein complexes are observed in solution, enabling the recording of images of the protein in multiple conformations. Various methods exist for capturing the conformational variability through analysis of cryo-EM data. Here, we analyzed the conformational variability in the hexameric AAA + ATPase p97, a complex with a six-fold rotational symmetric core surrounded by six flexible N-domains. We compared the performance of discrete classification methods with our recently developed method, MDSPACE, which uses 3D-to-2D flexible fitting of an atomic structure to images based on molecular dynamics (MD) simulations. Our analysis detected a novel conformation adopted by approximately 2% of the particles in the dataset and determined that the N-domains of p97 sway by up to 60° around a central position. This study demonstrates the application of MDSPACE in analyzing the continuous conformational changes in partially symmetrical protein complexes, systems notoriously difficult to analyze due to the alignment errors caused by their partial symmetry.

Structural insight into an Arl1-ArfGEF complex involved in Golgi recruitment of a GRIP-domain golgin.

Arl1 is an Arf-like (Arl) GTP-binding protein that interacts with the guanine nucleotide exchange factor Gea2 to recruit the golgin Imh1 to the Golgi. The Arl1-Gea2 complex also binds and activates the phosphatidylserine flippase Drs2 and these functions may be related, although the underlying molecular mechanism is unclear. Here we report high-resolution cryo-EM structures of the full-length Gea2 and the Arl1-Gea2 complex. Gea2 is a large protein with 1459 residues and is composed of six domains (DCB, HUS, SEC7, HDS1-3). We show that Gea2 assembles a stable dimer via an extensive interface involving hydrophobic and electrostatic interactions in the DCB and HUS region. Contrary to the previous report on a Gea2 homolog in which Arl1 binds to the dimerization surface of the DCB domain, implying a disrupted dimer upon Arl1 binding, we find that Arl1 binds to the outside surface of the Gea2 DCB domain, leaving the Gea2 dimer intact. The interaction between Arl1 and Gea2 involves the classic FWY aromatic residue triad as well as two Arl1-specific residues. We show that key mutations that disrupt the Arl1-Gea2 interaction abrogate Imh1 Golgi association. This work clarifies the Arl1-Gea2 interaction and improves our understanding of molecular events in the membrane trafficking.

Komagataella phaffii Erp41 is a protein disulfide isomerase with unprecedented disulfide bond catalyzing activity when coupled to glutathione.

In the methylotrophic yeast Komagataella phaffii, we identified an endoplasmic reticulum-resident protein disulfide isomerase (PDI) family member, Erp41, with a peculiar combination of active site motifs. Like fungal ERp38, it has two thioredoxin-like domains which contain active site motifs (a and a'), followed by an alpha-helical ERp29c C-terminal domain (c domain). However, while the a domain has a typical PDI-like active site motif (CGHC), the a' domain instead has CGYC, a glutaredoxin-like motif which confers to the protein an exceptional affinity for GSH/GSSG. This combination of active site motifs has so far been unreported in PDI-family members. Homology searches revealed ERp41 is present in the genome of some plants, fungal parasites, and a few nonconventional yeasts, among which are Komagataella spp. and Yarrowia lipolytica. These yeasts are both used for the production of secreted recombinant proteins. Here, we analyzed the activity of K. phaffii Erp41. We report that it is nonessential in K. phaffii, and that it can catalyze disulfide bond formation in partnership with the sulfhydryl oxidase Ero1 in vitro with higher turnover rates than the canonical PDI from K. phaffii, Pdi1, but slower activation times. We show how Erp41 has unusually fast glutathione-coupled oxidation activity and relate it to its unusual combination of active sites in its thioredoxin-like domains. We further describe how this determines its unusually efficient catalysis of dithiol oxidation in peptide and protein substrates.

Mechanism of Membrane Curvature Induced by SNX1: Insights from Molecular Dynamics Simulations.

SNX proteins have been found to induce membrane remodeling to facilitate the generation of transport carriers in endosomal pathways. However, the molecular mechanism of membrane bending and the role of lipids in the bending process remain elusive. Here, we conducted coarse-grained molecular dynamics simulations to investigate the role of the three structural modules (PX, BAR, and AH) of SNX1 and the PI3P lipids in membrane deformation. We observed that the presence of all three domains is essential for SNX1 to achieve a stable membrane deformation. BAR is capable of remodeling the membrane through the charged residues on its concave surface, but it requires PX and AH to establish stable membrane binding. AH penetrates into the lipid membrane, thereby promoting the induction of membrane curvature; however, it is inadequate on its own to maintain membrane bending. PI3P lipids are also indispensable for membrane remodeling, as they play a dominant role in the interactions of lipids with the BAR domain. Our results enhance the comprehension of the molecular mechanism underlying SNX1-induced membrane curvature and help future studies of curvature-inducing proteins.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes.

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 µM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

An unusual dual sugar-binding lectin domain controls the substrate specificity of a mucin-type O-glycosyltransferase.

N-acetylgalactosaminyl-transferases (GalNAc-Ts) initiate mucin-type O-glycosylation, an abundant and complex posttranslational modification that regulates host-microbe interactions, tissue development, and metabolism. GalNAc-Ts contain a lectin domain consisting of three homologous repeats (α, ß, and γ), where α and ß can potentially interact with O-GalNAc on substrates to enhance activity toward a nearby acceptor Thr/Ser. The ubiquitous isoenzyme GalNAc-T1 modulates heart development, immunity, and SARS-CoV-2 infectivity, but its substrates are largely unknown. Here, we show that both α and ß in GalNAc-T1 uniquely orchestrate the O-glycosylation of various glycopeptide substrates. The α repeat directs O-glycosylation to acceptor sites carboxyl-terminal to an existing GalNAc, while the ß repeat directs O-glycosylation to amino-terminal sites. In addition, GalNAc-T1 incorporates α and ß into various substrate binding modes to cooperatively increase the specificity toward an acceptor site located between two existing O-glycans. Our studies highlight a unique mechanism by which dual lectin repeats expand substrate specificity and provide crucial information for identifying the biological substrates of GalNAc-T1.

Structural analysis of PTPN21 reveals a dominant-negative effect of the FERM domain on its phosphatase activity.

PTPN21 belongs to the four-point-one, ezrin, radixin, moesin (FERM) domain-containing protein tyrosine phosphatases (PTP) and plays important roles in cytoskeleton-associated cellular processes like cell adhesion, motility, and cargo transport. Because of the presence of a WPE loop instead of a WPD loop in the phosphatase domain, it is often considered to lack phosphatase activity. However, many of PTPN21's biological functions require its catalytic activity. To reconcile these findings, we have determined the structures of individual PTPN21 FERM, PTP domains, and a complex between FERM-PTP. Combined with biochemical analysis, we have found that PTPN21 PTP is weakly active and is autoinhibited by association with its FERM domain. Disruption of FERM-PTP interaction results in enhanced ERK activation. The oncogenic HPV18 E7 protein binds to PTP at the same location as PTPN21 FERM, indicating that it may act by displacing the FERM domain from PTP. Our results provide mechanistic insight into PTPN21 and benefit functional studies of PTPN21-mediated processes.

Charged Amino Acids in the Transmembrane Helix Strongly Affect the Enzyme Activity of Aromatase.

Estrogens play critical roles in embryonic development, gonadal sex differentiation, behavior, and reproduction in vertebrates and in several human cancers. Estrogens are synthesized from testosterone and androstenedione by the endoplasmic reticulum membrane-bound P450 aromatase/cytochrome P450 oxidoreductase complex (CYP19/CPR). Here, we report the characterization of novel mammalian CYP19 isoforms encoded by CYP19 gene copies. These CYP19 isoforms are all defined by a combination of mutations in the N-terminal transmembrane helix (E42K, D43N) and in helix C of the catalytic domain (P146T, F147Y). The mutant CYP19 isoforms show increased androgen conversion due to the KN transmembrane helix. In addition, the TY substitutions in helix C result in a substrate preference for androstenedione. Our structural models suggest that CYP19 mutants may interact differently with the membrane (affecting substrate uptake) and with CPR (affecting electron transfer), providing structural clues for the catalytic differences.

Interaction of human dendritic cell receptor DEC205/CD205 with keratins.

DEC205 (CD205) is one of the major endocytic receptors on dendritic cells and has been widely used as a receptor target in immune therapies. It has been shown that DEC205 can recognize dead cells through keratins in a pH-dependent manner. However, the mechanism underlying the interaction between DEC205 and keratins remains unclear. Here we determine the crystal structures of an N-terminal fragment of human DEC205 (CysR∼CTLD3). The structural data show that DEC205 shares similar overall features with the other mannose receptor family members such as the mannose receptor and Endo180, but the individual domains of DEC205 in the crystal structure exhibit distinct structural features that may lead to specific ligand binding properties of the molecule. Among them, CTLD3 of DEC205 adopts a unique fold of CTLD, which may correlate with the binding of keratins. Furthermore, we examine the interaction of DEC205 with keratins by mutagenesis and biochemical assays based on the structural information and identify an XGGGX motif on keratins that can be recognized by DEC205, thereby providing insights into the interaction between DEC205 and keratins. Overall, these findings not only improve the understanding of the diverse ligand specificities of the mannose receptor family members at the molecular level but may also give clues for the interactions of keratins with their binding partners in the corresponding pathways.

Identification of a conserved α-helical domain at the N terminus of human DNA methyltransferase 1.

In vertebrates, DNA methyltransferase 1 (DNMT1) contributes to preserving DNA methylation patterns, ensuring the stability and heritability of epigenetic marks important for gene expression regulation and the maintenance of cellular identity. Previous structural studies have elucidated the catalytic mechanism of DNMT1 and its specific recognition of hemimethylated DNA. Here, using solution nuclear magnetic resonance spectroscopy and small-angle X-ray scattering, we demonstrate that the N-terminal region of human DNMT1, while flexible, encompasses a conserved globular domain with a novel α-helical bundle-like fold. This work expands our understanding of the structure and dynamics of DNMT1 and provides a structural framework for future functional studies in relation with this new domain.

Deciphering the stereo-specific catalytic mechanisms of cis-epoxysuccinate hydrolases producing L(+)-tartaric acid.

Microbial epoxide hydrolases, cis-epoxysuccinate hydrolases (CESHs), have been utilized for commercial production of enantiomerically pure L(+)- and D(-)-tartaric acids for decades. However, the stereo-catalytic mechanism of CESH producing L(+)-tartaric acid (CESH[L]) remains unclear. Herein, the crystal structures of two CESH[L]s in ligand-free, product-complexed, and catalytic intermediate forms were determined. These structures revealed the unique specific binding mode for the mirror-symmetric substrate, an active catalytic triad consisting of Asp-His-Glu, and an arginine providing a proton to the oxirane oxygen to facilitate the epoxide ring-opening reaction, which has been pursued for decades. These results provide the structural basis for the rational engineering of these industrial biocatalysts.

FERM domain-containing proteins are active components of the cell nucleus.

The FERM domain is a conserved and widespread protein module that appeared in the common ancestor of amoebae, fungi, and animals, and is therefore now found in a wide variety of species. The primary function of the FERM domain is localizing to the plasma membrane through binding lipids and proteins of the membrane; thus, for a long time, FERM domain-containing proteins (FDCPs) were considered exclusively cytoskeletal. Although their role in the cytoplasm has been extensively studied, the recent discovery of the presence and importance of cytoskeletal proteins in the nucleus suggests that FDCPs might also play an important role in nuclear function. In this review, we collected data on their nuclear localization, transport, and possible functions, which are still scattered throughout the literature, with special regard to the role of the FERM domain in these processes. With this, we would like to draw attention to the exciting, new dimension of the role of FDCPs, their nuclear activity, which could be an interesting novel direction for future research.

STK19 is a DNA/RNA-binding protein critical for DNA damage repair and cell proliferation.

STK19 was originally identified as a manganese-dependent serine/threonine-specific protein kinase, but its function has been highly debated. Here, the crystal structure of STK19 revealed that it does not contain a kinase domain, but three intimately packed winged helix (WH) domains. The third WH domain mediated homodimerization and double-stranded DNA binding, both being important for its nuclear localization. STK19 participated in the nucleotide excision repair (NER) and mismatch repair (MMR) pathways by recruiting damage repair factors such as RPA2 and PCNA. STK19 also bound double-stranded RNA through the DNA-binding interface and regulated the expression levels of many mRNAs. Furthermore, STK19 knockdown cells exhibited very slow cell proliferation, which cannot be rescued by dimerization or DNA-binding mutants. Therefore, this work concludes that STK19 is highly unlikely to be a kinase but a DNA/RNA-binding protein critical for DNA damage repair (DDR) and cell proliferation. To prevent further confusions, we renamed this protein as TWH19 (Tandem Winged Helix protein formerly known as STK19).

Methylation of elongation factor 1A by yeast Efm4 or human eEF1A-KMT2 involves a beta-hairpin recognition motif and crosstalks with phosphorylation.

Translation elongation factor 1A (eEF1A) is an essential and highly conserved protein required for protein synthesis in eukaryotes. In both Saccharomyces cerevisiae and human, five different methyltransferases methylate specific residues on eEF1A, making eEF1A the eukaryotic protein targeted by the highest number of dedicated methyltransferases after histone H3. eEF1A methyltransferases are highly selective enzymes, only targeting eEF1A and each targeting just one or two specific residues in eEF1A. However, the mechanism of this selectivity remains poorly understood. To reveal how S. cerevisiae elongation factor methyltransferase 4 (Efm4) specifically methylates eEF1A at K316, we have used AlphaFold-Multimer modeling in combination with crosslinking mass spectrometry (XL-MS) and enzyme mutagenesis. We find that a unique beta-hairpin motif, which extends out from the core methyltransferase fold, is important for the methylation of eEF1A K316 in vitro. An alanine mutation of a single residue on this beta-hairpin, F212, significantly reduces Efm4 activity in vitro and in yeast cells. We show that the equivalent residue in human eEF1A-KMT2 (METTL10), F220, is also important for its activity towards eEF1A in vitro. We further show that the eEF1A guanine nucleotide exchange factor, eEF1Bα, inhibits Efm4 methylation of eEF1A in vitro, likely due to competitive binding. Lastly, we find that phosphorylation of eEF1A at S314 negatively crosstalks with Efm4-mediated methylation of K316. Our findings demonstrate how protein methyltransferases can be highly selective towards a single residue on a single protein in the cell.

An imine reductase that captures reactive intermediates in the biosynthesis of the indolocarbazole reductasporine.

Imine reductases (IREDs) and reductive aminases have been used in the synthesis of chiral amine products for drug manufacturing; however, little is known about their biological contexts. Here we employ structural studies and site-directed mutagenesis to interrogate the mechanism of the IRED RedE from the biosynthetic pathway to the indolocarbazole natural product reductasporine. Cocrystal structures with the substrate-mimic arcyriaflavin A reveal an extended active site cleft capable of binding two indolocarbazole molecules. Site-directed mutagenesis of a conserved aspartate in the primary binding site reveals a new role for this residue in anchoring the substrate above the NADPH cofactor. Variants targeting the secondary binding site greatly reduce catalytic efficiency, while accumulating oxidized side-products. As indolocarbazole biosynthetic intermediates are susceptible to spontaneous oxidation, we propose the secondary site acts to protect against autooxidation, and the primary site drives catalysis through precise substrate orientation and desolvation effects. The structure of RedE with its extended active site can be the starting point as a new scaffold for engineering IREDs and reductive aminases to intercept large substrates relevant to industrial applications.

Subdomain dynamics enable chemical chain reactions in non-ribosomal peptide synthetases.

Many peptide-derived natural products are produced by non-ribosomal peptide synthetases (NRPSs) in an assembly-line fashion. Each amino acid is coupled to a designated peptidyl carrier protein (PCP) through two distinct reactions catalysed sequentially by the single active site of the adenylation domain (A-domain). Accumulating evidence suggests that large-amplitude structural changes occur in different NRPS states; yet how these molecular machines orchestrate such biochemical sequences has remained elusive. Here, using single-molecule Förster resonance energy transfer, we show that the A-domain of gramicidin S synthetase I adopts structurally extended and functionally obligatory conformations for alternating between adenylation and thioester-formation structures during enzymatic cycles. Complementary biochemical, computational and small-angle X-ray scattering studies reveal interconversion among these three conformations as intrinsic and hierarchical where intra-A-domain organizations propagate to remodel inter-A-PCP didomain configurations during catalysis. The tight kinetic coupling between structural transitions and enzymatic transformations is quantified, and how the gramicidin S synthetase I A-domain utilizes its inherent conformational dynamics to drive directional biosynthesis with a flexibly linked PCP domain is revealed.

A conserved electrostatic membrane-binding surface in synaptotagmin-like proteins revealed using molecular phylogenetic analysis and homology modeling.

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among the vertebrate synaptotagmin-like protein (Slp) family. Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

Klebsiella pneumoniae carbapenemase variant 44 acquires ceftazidime-avibactam resistance by altering the conformation of active-site loops.

Klebsiella pneumoniae carbapenemase 2 (KPC-2) is an important source of drug resistance as it can hydrolyze and inactivate virtually all ß-lactam antibiotics. KPC-2 is potently inhibited by avibactam via formation of a reversible carbamyl linkage of the inhibitor with the catalytic serine of the enzyme. However, the use of avibactam in combination with ceftazidime (CAZ-AVI) has led to the emergence of CAZ-AVI-resistant variants of KPC-2 in clinical settings. One such variant, KPC-44, bears a 15 amino acid duplication in one of the active-site loops (270-loop). Here, we show that the KPC-44 variant exhibits higher catalytic efficiency in hydrolyzing ceftazidime, lower efficiency toward imipenem and meropenem, and a similar efficiency in hydrolyzing ampicillin, than the WT KPC-2 enzyme. In addition, the KPC-44 variant enzyme exhibits 12-fold lower AVI carbamylation efficiency than the KPC-2 enzyme. An X-ray crystal structure of KPC-44 showed that the 15 amino acid duplication results in an extended and partially disordered 270-loop and also changes the conformation of the adjacent 240-loop, which in turn has altered interactions with the active-site omega loop. Furthermore, a structure of KPC-44 with avibactam revealed that formation of the covalent complex results in further disorder in the 270-loop, suggesting that rearrangement of the 270-loop of KPC-44 facilitates AVI carbamylation. These results suggest that the duplication of 15 amino acids in the KPC-44 enzyme leads to resistance to CAZ-AVI by modulating the stability and conformation of the 270-, 240-, and omega-loops.

Structure and biochemical characterization of l-2-hydroxyglutarate dehydrogenase and its role in the pathogenesis of l-2-hydroxyglutaric aciduria.

l-2-hydroxyglutarate dehydrogenase (L2HGDH) is a mitochondrial membrane-associated metabolic enzyme, which catalyzes the oxidation of l-2-hydroxyglutarate (l-2-HG) to 2-oxoglutarate (2-OG). Mutations in human L2HGDH lead to abnormal accumulation of l-2-HG, which causes a neurometabolic disorder named l-2-hydroxyglutaric aciduria (l-2-HGA). Here, we report the crystal structures of Drosophila melanogaster L2HGDH (dmL2HGDH) in FAD-bound form and in complex with FAD and 2-OG and show that dmL2HGDH exhibits high activity and substrate specificity for l-2-HG. dmL2HGDH consists of an FAD-binding domain and a substrate-binding domain, and the active site is located at the interface of the two domains with 2-OG binding to the re-face of the isoalloxazine moiety of FAD. Mutagenesis and activity assay confirmed the functional roles of key residues involved in the substrate binding and catalytic reaction and showed that most of the mutations of dmL2HGDH equivalent to l-2-HGA-associated mutations of human L2HGDH led to complete loss of the activity. The structural and biochemical data together reveal the molecular basis for the substrate specificity and catalytic mechanism of L2HGDH and provide insights into the functional roles of human L2HGDH mutations in the pathogeneses of l-2-HGA.